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THE PERTURBATIVE APPROACH TO PATH INTEGRALS: 
A SUCCINCT MATHEMATICAL TREATMENT 

TIMOTHY NGUYEN 


Abstract. We study finite-dimensional integrals in a way that elucidates the mathemat¬ 
ical meaning behind the formal manipulations of path integrals occurring in quantum field 
theory. This involves a proper understanding of how Wick’s theorem allows one to evaluate 
integrals perturbatively, i.e., as a series expansion in a formal parameter irrespective of 
convergence properties. We establish invariance properties of such a Wick expansion under 
coordinate changes and the action of a Lie group of symmetries, and we use this to study 
essential features of path integral manipulations, including coordinate changes. Ward iden¬ 
tities, Schwinger-Dyson equations, Faddeev-Popov gauge-fixing, and eliminating fields by 
their equation of motion. We also discuss the asymptotic nature of the Wick expansion and 
the implications this has for defining path integrals perturbatively and nonperturbatively. 
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Introduction 

Quantum field theory often makes use of manipulations of path integrals that are without 
a proper mathematical definition and hence have only a formal meaning. These computa¬ 
tions are usually presented in a manner that makes it difficult to ascertain the mathematical 
nature of these operations, other than that they are inspired from familiar properties of 
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finite-dimensional integrals. The purpose of this paper is to explain in a large class of im¬ 
portant examples how formal path integral manipulations constitute notation for otherwise 
mathematically well-defined procedures. While some practitioners of quantum field theory 
are aware of these issues (see e.g. [23]), considerations of rigor are highly non-uniform in 
the quantum field theoretic literature, making it a challenge to navigate for both beginner 
and expert alike. The present note arose from the need to have a simple explanation of 
the mathematics behind perturbative quantum field theory, one that is both rigorous and 
sufficiently comprehensive to cover examples of interest. 

A basic prototype for the kinds of integrals one studies in quantum mechanics and quan¬ 
tum field theory is the one-dimensional integral 


I{X) = 


dx e 


Ax 


n 


( 0 . 1 ) 


with A a real parameter and n > 3. The standard procedure of expanding e as a Taylor 

1 2 

series and then integrating term by term against the Gaussian measure dxe~ 2 ^ yields a 
well-defined formal power series in A irrespective of whether /(A) converges (which holds 
only for A > 0 when n is even and A = 0 when n is odd). In particular, we can formally 
manipulate the integrand occurring in /(A), such as performing a change of variables, and 
still produce a corresponding formal series in A using the same method. 

This procedure of evaluating an integral as a formal power series in the relevant coupling 
constants is the guiding principle behind perturbative quantum field theory. In this paper, 
we refer to such a procedure as the Wick expansion. It arises from a formal application of 
the saddle point approximation and does not require the underlying integral to converge. In 
the infinite-dimensional setting of quantum field theory, the Wick expansion appropriately 
generalized (which leads to the familiar regularization and renormalization of Feynman 
diagrams) provides a definition of path integrals. Such a definition bypasses the difficulties 
(or even impossibility) of constructing the appropriate measures on infinite-dimensional 
spaces needed to evaluate a path integral as an honest integral. Nevertheless, one refers 
to the path integral as being defined perturbatively through this formal series method, as 
though the series were known to approximate path integrals as honest integrals. (This is 
inspired from the fact that for n even, the Wick expansion for /(A) yields an asymptotic 
series for /(A), regarded as an honest function of A, as A —)• O"''.) As a consequence, the 
many results in quantum field theory obtained from formal manipulations of path integrals 
naturally arouse suspicion. However, based on the finite-dimensional case above, the way 
out of this confusion is to have a clear separation between integration, which is analytic, 
and the Wick expansion, which is algebraic. 

This paper provides a mathematical study of the Wick expansion of finite-dimensional 
integrals in a way that provides insight into (perturbative) path integrals in quantum field 
theory. Namely, we establish invariance properties of Wick expansions with respect to 
standard calculus manipulations of the integrals used to define them and show how this 
allows one to make mathematical meaning of formal path integral manipulations in quantum 
field theory. It is helpful to illustrate this approach via the following schema: 
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Finite-Dimensional Integrals 


1(A) 


Wick expansion 


formal series 
in A 

-r 


integration 

¥ 

function of A 


asymptotic series 


Perturbative Path Integrals 




Wick expansion 


formal series 
in A 


Figure 1. The Wick expansion and its relationship to integrals. 


Thus, while integral manipulations can be analyzed in the usual way when integration is 
well-defined, they can also (and always) be analyzed using the Wick expansion. 

This paper is organized as follows. In Section 1, we define the Wick expansion and 
show that it is independent of the choice of coordinates used to define it. The latter is 
important for showing that the Wick expansion can be defined on manifolds. It should 
be emphasized that the Wick expansion is an expansion about an arbitrary nondegenerate 
critical point (the Hessian can be indefinite), so that such an expansion is not always a 
Gaussian approximation. We also study the more general Morse-Bott case in which we 
Wick expand about critical submanifolds. In Section 2, we discuss the Wick expansion in 
the presence of symmetries. Here, a suitable gauge-fixing procedure is needed to define 
the Wick expansion. In doing so, we obtain a rigorous formulation of the Faddeev-Popov 
gauge-fixing procedure in the finite-dimensional setting, which works whenever we have a 
volume preserving action of a (not necessarily compact) Lie group on a manifold. In fact, we 
emphasize the there are two types of Faddeev-Popov procedures (we call them the slice and 
weighted versions), and we clarify both of these in the context of integration and the Wick 
expansion in ways that are crucial for understanding quantum gauge theories. In Section 
3, we relate the Wick expansion to the asymptotics of integrals (this is the diagonal arrow 
in Figure 1), thereby relating purely formal algebraic manipulations to analytic properties 
of integrals which are convergent and satisfy additional hypotheses. This provides the 
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justification for why the Wick expansion is a sensible algebraic object to associate to an 
integral. Finally in Section 4 we discuss how to interpret our results in the context of 
quantum field theory. Namely, we discuss the Wick expansion in light of the regularization 
procedures typically used in quantum field theory - either from regulating the propagator, 
the dimension, or using a lattice. From this, we can suitably interpret a variety of formal 
manipulations of path integrals defined perturbatively in a mathematically rigorous manner. 
By comparison, standard textbook treatments of these procedures (see e.g. [m Ch 16.4] 
|20l Ch 15.7]) interpret them via integration, which is often illegitimate (both literally in 
the infinite-dimensional setting and by analogy with the finite-dimensional setting). This is 
particularly true of gauge-fixing procedures, for which the blur between integration and the 
Wick expansion hampers a rigorous understanding as our discussion will show. We believe 
the examples we discuss are only a few among many that achieve clarity (and perhaps even 
conceptual correction) through the methodology presented in this paper, which is centered 
around the Wick expansion. We conclude with a brief discussion concerning perturbative 
versus constructive quantum field theory. 

Altogether, our treatment separates in a simple and elegant manner the analytic difficul¬ 
ties of working in infinite dimensions from the distinction between integration and the Wick 
expansion already occurring in finite dimensions. In other words, we clearly distinguish be¬ 
tween analysis (integration) and algebra (Wick expansion). Typically, analytic approaches 
to quantum field theory (e.g. those arising from Wightman and Osterwalder-Schrader ax¬ 
ioms HH) do not resemble the treatment of perturbative quantum field theory that appears 
in textbooks [S] fTHl 120] . whose methods have proven to be remarkably successfnl in real 
world applications. On the other hand, algebraic approaches to quantum field theory are 
often limited to or are only concerned with producing formal series expansions, which even 
in finite dimensions, do not capture integration withont further analytical considerations. 
We hope our explicit delineation of algebraic and analytic methods provides clarity to the 
limitations that various approaches to quantum field theory have. Perhaps a helpful analogy 
would be to regard algebraic methods that produce a formal perturbative series as providing 
a “weak” construction of a path integral; analytic methods which attempt to boost such a 
formal series to a well-defined function then provide a “strong” construction. As is true in 
many other contexts, it is conceptually sound to have a clear distinction between weak and 
strong constructions of an object. 

In order to keep this note succint, we wrote it as a compromise between comprehensiveness 
and brevity with emphasis on simplicity. For alternative treatments to some of the issues 
raised here, see e.g. [laiiT!. 


1. The Wick Expansion 


Consider the integral 

1 = J 


( 1 . 1 ) 


on where S and / are complex-valued functions. We always assume that our functions 
are smooth, i.e., infinitely differentiable. For the situations relevant to physics, we regard 
S{x) as an action, f{x) an observable, and the parameter h a real or complex parameter. 
Moreover, we think of h as being small, since li —>■ 0 is to be regarded as a semiclassical 
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limit of (1.1). A rescaling of (1.1) shows that h can also be regarded as a perturbative 


parameter, since for example, letting S{x) = + x^, f{x) = 1, and x i—)• we recover 

the integral considered in the introduction with A = 

Under the appropriate hypotheses, the integral I is convergent and becomes a function 
of h. On the other hand, one can treat h as a purely formal parameter, in which case I is 
a formal object consisting of an integral sign and an integrand. One can then attempt to 
dehne various series expansions in h from I using purely algebraic rules. One such expansion 
is the Wick expansion, which takes as an additional source of input a nondegenerate critical 
point xq of S{x). Under appropriate additional assumptions (see Section 3), the Wick 
expansion about xq provides the asymptotics of I when /(x) is a bump function localized 
around a nondegenerate critical point of S{x). However, the Wick expansion is always a 
well-defined series regardless of the convergence properties of I, since it depends only on 
the derivatives of the integrand of I at xq. 

In fact, the significance of the algebraic nature of the Wick expansion goes beyond just 
bypassing convergence issues of /, since the latter can always be made convergent by making 
/(x) compactly supported. In the infinite dimensional setting of quantum held theory, an 
honest integral I is replaced with a formal path integral. One often does not know how 
to make sense of the integrand of such a path integral, since it may lack a well-dehned 
construction as an honest measure. However, the Wick expansion, suitably interpreted, 
provides an algebraic way to integrate formally such an ill-dehned measure. This is based 
on the analogy between the Wick expansion and integration that will become apparent in 
what follows. We return to the bearing these observations have on quantum held theory in 
Section 4. 

We now proceed to dehne the Wick expansion. Consider a nondegenerate critical point 
Xq of S'(x), i.e., one for which the Hessian of S{x) at xq is nondegenerate. Explicitly, we 
can write 


1 

S{x) = 5(xo) -A{x - xo, X - xo) -h 0(|x - xo| ) 


( 1 . 2 ) 


for X near xq, where the symmetric bilinear pairing A = A(-, •) is the Hessian. Thus, we 
can write 


^-S(x)/h _ ^-S(xo)/h^-A{x-xo,x-xo)/2h 


oS(x)/h 


(1.3) 


the product of a Gaussian and an “interaction” term. (For the time being, we assume 
A is positive dehnite, though we will remove this assumption shortly.) In (1.1), we can 
expand as a Taylor series centered at xq and then integrate term by term against 

the Gaussian measure e~^^^AI^Q-M^-^o,x-xo))/ 2 h^ Aside from an overall normalization 

factor, we obtain a formal power series in h. Indeed, one can easily see this from the 
rescaling x i-A- xq + 1t}^‘^{x — xq). This makes the Gaussian measure proportional to 
and then the integration of polynomials of even degree only picks up even powers of 

The above construction yields the Wick expansion of I about xq in case A is positive 
dehnite. However, it can be generalized to A nondegenerate as follows. The hrst step is 
to recall a result that goes by the name of Wick’s Theorem, which provides a convenient, 
combinatorial formula for evaluating integrals of polynomials against Gaussian measures. 
The easiest way to describe a polynomial function on a vector space is to pick a basis. On 
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we can pick a basis ei,..., and a corresponding dual basis • • • , of coordinate 
monomials. We abbreviate 

d'^x = dx^ ■ ■ ■ dx^ 

for the corresponding density, which we assume to be the same density that appears in 
We can describe a normalized Gaussian measure in terms of the matrix 


via 


dfj,A = dx^ ■ ■ ■ dx'^ ^ 


det^ij^ ^ ^-A{x,x)l2 


(271) 


Let denote the inverse matrix of Aij. 


Theorem 1.1. (Wick’s Theorem) For A positive definite, we have 


dp,A X* 




2™m! ^ 

0'€iS2rr 


<T(2m —l)^cr(2m) 


(1.4) 


This formula can be encoded pictorially through the use of Feynman diagrams, for which 
the contractions of the A’'^ into the slots of x*^ • • • x*^™ are encoded through incidence re¬ 
lations among edges and vertices of graphs. Further details can be found e.g. in UM or 
other textbooks on quantum field theory. 


Note that while the left-hand side of (1.4) is analytic in nature, the right-hand side is 


purely combinatorial. In particular, one can define the right-hand side for any nonsingular 


(complex-valued) matrix A. One can also regard the right-hand side of (1.4) as providing 


an analytic continuation of the left-hand side to the space of nonsingular matrices. Hence, 
we make the following definition: 


Definition 1.2. Given a nondegenerate, symmetric bilinear pairing A, define the Wick 
operator Wa to be the linear functional on the space of polynomials given by the following 
formula; 


Wa(P) 


1 


2™m! 


y ylMi)V(2)... ^* a (2m—1)^17 (2m) 


crE5'2r) 


0 


P is odd. 


(1.5) 


Although the above formula makes use of a basis, it is easy to see that the definition of 
WA depends only on the bilinear pairing A. The terms appearing in the Wick formula 
(1.5) are called Wick contractions. 

Wick’s Theorem tells us that the integration of polynomials against the Gaussian measure 
dpA coincides with the Wick operator Wa for A positive definite. Thus, the series expansion 
which we described above, which involved integrating the Taylor series of /(x)e'^^’^^/^ against 
the Gaussian can be defined using the Wick operator Wa 

instead of integration. But since the Wick operator is well-defined for any nondegenerate 
A, this allows us to extend the definition of such an expansion to the case when S{x) has 
a critical point with arbitrary nondegenerate Hessian; 
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Definition 1.3. Consider the integral I along with a choice of local coordinates near a 
nondegenerate critical point xq of S{x). This yields a fixed decomposition of the integrand 
of I into the product of t he c oordinate density df^x and the function near xq. 

Represent as in (1.3) and Taylor expand 

/(xo + h^/‘^{x - 

about X = xo to obtain a power series Yl^=o “ ®o) grouped by powers of h. Then 

the Wick expansion of I about xq is the formal series in h given by 

1/2 


IT. 


“ ' det A 


{2TrhY 




^-S{xo)/h 




(1.7) 


k=0 


where Ck = WAiPk)- 


A priori, the Wick expansion depends on the choice of coordinates used to define it. 
Indeed, a different choice of coordinates leads to a Jacobian factor and a different sequence 


of polynomials occurring in the Taylor expansion of (1.6). Nevertheless, Theorem 1.5 tells 


us that the Wick expansion yields a series independent of the coordinate system used to 
construct it. First, we prove the following lemma: 

Lemma 1.4. The Wick expansion of a total derivative is zero. 

Proof We need to show that the Wick expansion of 




= K 


-1 


d'^xe ^^^">/\hd^if{x) - d^iS{x)f{x)] 


( 1 . 8 ) 


vanishes. Here, the Wick expansion is obtained by applying Definition 1.3 to t he f inal 
expression above (the h~^ is just an overall constant factor). Using the splitting (1.3), it 


follows that in the case of A positive definite, the vanishing of the Wick expansion follows 
from 

J dh^xd^r [p(x)e-^(^’^)/2j (19) 

for arbitrary polynomials P{x). For general nondegenerate A, one has to establish the 
algebraic analogue of the above equation, namely 

Wa {dxiP{x) — AijXjP{x)) = 0. (1-10) 

Verifying this identity is straightforward. □ 


Theorem 1.5. The Wick expansion about a nondegenerate critical point xq is independent 
of the choice of coordinates. 

Proof Let —>■ be a diffeomorphism. We want to show that the Wick expansion 

of d^^x f{x)e~^^P/^ and {d'^x f{x)e~^P'>/^) about xq and <I>“^(xo), respectively, are equal 
as series expansions in h. Without loss of generality, we can suppose xq = 0 and $(0) = 0. 
Moreover, it is easy to see that the Wick expansion is invariant under a linear change 
of coordinates, since linear maps preserve polynomial degree and hence the Wick formula 
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(1.5). Thus, we may further suppose that at x = 0 is the identity. It follows that the 
1-parameter family of maps 

<l>t(x) = (1 — t)x -)- t^{x) 

is a local diffeomorphism in a neighborhood of 0 for all t € [0,1]. It suffices to show that 
the Wick expansion of ^^'Kd^x f{x)e~^^^^^^) is identically zero for all t. Letting 


= is 


(x)) 


( 1 . 11 ) 


S=t 


denote the time-dependent vector field associated to the flow then 






= di, 




<^>1 (d“'x/(x)e-^(^)/^ 


( 1 . 12 ) 


where Ty is the Lie derivative, and in the second line we used the Cartan formula £y = dty 
for the Lie derivative of a differential form of top degree. Since <I>t(x) = x -|- 0(|xp) for all 
t, then the i$*(y)<l>| {d'^x f{x)e~^^^^/^') are all of the form for 

some t-dependent functions St and ft- Hence, the integrands (1.12) are all integrands for 
which we may Wick expand. By Lemma 1.4, these Wick expansions are all zero. □ 


Remark 1.6. The Wick expansion only depends on the derivatives of the integrand of I 
at xq, i.e., the infinite jet. Thus, one should really work in the category of formal power 
series about xq (e.g. changes of coordinates need only be invertible formal power series). 
We will leave this understanding implicit and instead maintain more geometric terminology 
throughout the paper by considering all objects as smooth. 


1.1. The Morse-Bott case. The coordinate invariance of the Wick expansion allows us 
to generalize our definition of it on Euclidean space to the setting of smooth manifolds. 
Consider the integral 

1= [ dVf{x)e-^^^'>/^ (1.13) 

Jm 

where dV is a differential form of top degre^ on M. Note that the presence of /(x) is 
redundant since it can be grouped with dV, but we keep /(x) separate since one often keeps 
some reference dV fixed while varying /(x). Moreover, the presence of the h parameter in 
the exponential means that the (^-independent) function S'(x) is uniquely-defined, and we 
can refer to critical sets for ^(x) as critical sets of (the integrand of) /. Recall that a 
submanifold Z C M is critical for S if {dS)x = 0 for all x G Z. 

When Z = Xq is a nondegenerate critical point, we can work in local coordinates near 
Xq and Wick expand about xq. This yields a well-defined series in h independent of the 
coordinate system chosen by Theorem 1.5, We now want to consider the case when Z is 


^Our discussion generalizes straightforwardly to (and ought to be phrased in terms of) densities, but to 
keep our discussion a bit simpler and more familiar, we restrict ourselves to working with differential forms 
on orientable manifolds. 
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Morse-Bott nondegenerate, i.e., for every x £ Z, the kernel of the Hessian of S' at x equals 
TxZ. It is important to make this generalization since the types of integrals one often 
encounters are those for which S{x) has some moduli of stationary configurations. 

Given any Morse-Bott nondegenerate critical submanifold Z C M, consider a small 
tubular neighborhood Z oi Z. By regarding Z as a fibration over Z by disks (say by picking 
a Riemannian metric on M and using the exponential map in the direction orthogonal to Z), 
then the corresponding bundle projection vr ; Z —)• Z allows us to perform fiber integration 
[T]. This is a map tt* sending top-degree differential forms on Z to top-degree differential 
forms on Z, which is essentially integration along the fiber directions of Z. Moreover, fiber 
integration is volume-preserving: 



(1.14) 


Given the analogy between the Wick expansion and integration, we define a Wick expan¬ 
sion version of fiber integration, which we call Wick fiber integration. As one might expect, 
Wick fiber integration maps top-degree forms on Z to top-degree forms on Z valued in a 
formal series in h. This operation depends only on the local behavior near Z and hence 
extends to a map on top-degree forms on M. 

We define Wick fiber integration as follows. Given the bundle projection vr, we can 
choose local coordinates adapted to tt in the following sense. Given p G Z, we can choose 
coordinates {x,y) = {x^,... ,x^,y^,... on an open set U containing p such that 

ZnU = {y^ = ... = y'^~^ = 0} and vriZriG—>-Znf7is projection onto the x-coordinates. 
We refer to x and y as horizontal and vertical directions, respectively. The Morse-Bott 
nondegeneracy condition on Z implies that in these coordinates 


S{x, y) = S{x, 0) ]^Ax{y, y) + Sx{y), 


(1.15) 


where Ax{y, y) is a nondegenerate x-dependent quadratic form of the v ariabl es y^,..., y'^~'' 
and Sx{y) = 0{y^). Thus, with p £ Z fixed, we can Wick expand (1.13) about p with 
respect to the fiber variables y^,... ,y‘^~^, thereby obtaining a series in h valued in a top- 
degree differential form on Z. Indeed, the Wick expansion eliminates the y-variables and 
leaves the x-variables remaining. 


Abbreviate the integrand of (1.13) by 


V = dVf{x)e-^^^^/^. 


(1.16) 


Definition 1.7. Let Z C M be a submanifold and choose a fiber bundle structure vr : 
Z ^ Z. Wick fiber integration vr^ is the map which sends V, a top-degree differential form 
) with Z a Morse-Bott nondegenerate critical set, to (V), a top-degree differential 
form on Z valued in a formal series in h, by performing fiberwise Wick expansion of V along 
each point of Z. 



As before, we want to show that our definition of Wick fiber integration is independent 
of the choice of coordinates used to perform the Wick expansion along individual fibers. 

Lemma 1.8. The Wick fiber integration map is well-defined, i.e., it is independent of 
the choice of coordinates used to define it. 
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Proof Given coordinates {x, y) adapted to tt, another coordinate system adapted to tt 
is obtained via a local diffeomorphism y)), where and are the horizontal 

and vertical components of <h. Note that does not depend on y, since preserves the 
fibers of tt. The Wick expansion only depends on the Taylor series of all objects involved in 
y-variables. Thus, using Theorem 1.5, we can suppose ^'^{x,y) is the identity map on the 
y-variables up to translation, i.e. = ^'"{x) + y. In other words, maps the horizon¬ 

tal slice Z = {{x, 0)} to the set of points {(<h^(x), $^(x))} and it vertically translates fibers 
accordingly. Any such translation diffeomorphism commutes with the Wick expansion in 
the y-variables, i.e., one can Wick expand then push-forward by or else pushforward by 
and then Wick expand. This shows that vr^ is well-defined, since it transforms appro¬ 
priately under diffeomorphisms. □ 


Another way of reformulating the above is the following: 

Lemma 1.9. Wick fiber integration is covariant. That is, if ^ : M - 
morphism intertwining fiber bundle structures vr and <h*7r := <I) o vr o 
respectively, then 

:7rr(V)) = ($,7r)f(cI.,(V)). 




-> M is a diffeo- 
for Z and ^{Z), 

(1.17) 


Fiber integration maps top-degree differential forms on a total space to top-degree differ¬ 
ential forms on the base such that the resulting map on cohomology is independent of the 
homotopy class of the fiber bundle structure. This is because if = vr o , where is a 
one-parameter family of diffeomorphisms, then letting Vt denote the vector field associated 


to the flow as in (1.11), we have 




— TT* 


o <I)^ o diVt 
= d O TT* o <I> 


OLVf 


(1.18) 


Here, we used that fiber integration commutes with the exterior derivative. 

Wick fiber integration, being a formal analogue of fiber integration, satisfies the following 
property: 

Lemma 1.10. The terms of the Wick fiber integration vr^ (V) have well-defined cohomology 
classes that are independent of the choice of fiber bundle structure vr. 

Proof Given a smooth family vr^ of bundle structures, then by repeating the derivation 


(1.18) with Wick fiber integration in place of fiber integration, we have that (V) at 

p ^ Z consists of the Wick expansion along 7r^“^(p) of a total derivative. Near Z, d splits into 
a (time-dependent) fiber component and a component dz tangential to Z. Thus ^("t*)!^(V) 
is equal to a total dz derivative, since the fiber component of d is annihilated (by Lemma 
1.4) and dz commutes with Wick expansion in the fiber directions. This establishes the 


theorem for all fiber bundle structures homotopic to vr. Now we recall Remark 1.6 and note 


that Wick hber integration with respect to vr only depends on the Taylor expansion of vr 
in the directions transverse to Z (in particular, with respect to the fiber directions with 
respect to some fixed bundle structure); in the Morse-Bott situation, this Taylor expansion 
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has coefficients that are functions on Z. It is thus enough to show that the fiberwise Taylor 
expansions along Z of any two fiber bundle structures can be homotoped. 

To see this, note that given a fixed tt : Z —)• Z, we can always regard tt as a vector 
bundle projection tt : NZ ^ Z (where NZ is the normal bundle to Z with respect to 
some arbitrary metric on M) restricted to Z C NZ. With this normalization, the Taylor 
series of vr along the fibers of NZ is trivial. A general fiber bundle map vr' : Z —?■ Z has a 
nontrivial Taylor series along the fibers of NZ, since the latter fibers will not coincide with 
the fibers of tt' . The derivatives along Z of tt' in the fiber directions of NZ are sections of 
the bundle nA:>i This is a contractible space, with the zero section 

corresponding to vr. It follows that we can homotope the fiberwise Taylor expansions of any 
two fiber bundle structures. □ 


As formula (1.14) shows, fiber integration along Z splits integration in a neighborhood 
of Z into a fiber direction and a remaining direction parallel to Z. Our definition of the 
Wick expansion in the Morse-Bott setting is defined accordingly. Namely, the Wick expan¬ 
sion about Z involves a Wick fiber integration in the fiber directions, and then a residual 
integration along Z. More precisely, we have the following: 


Definition 1.11. Suppose Z is a compact Morse-Bott nondegenerate critical set for V. 
Define the Wick expansion of I about Z to be the formal series in h defined by 


Wz{n) = f Trf (dVf{x)e-^^^'^/^ 


(1.19) 


where tt is any fiber bundle structure for Z. Lemma 1.10| implies this definition is indepen¬ 
dent of the choice of vr. 


We emphasize that the integration over Z in the above is an honest integration, in which 
case we must make some assumptions to ensure that such an integration is well-defined. 
One cannot perform a residual Wick expansion on Z, since vr^ (^dV is of the 

form times a power series in h, with C the constant value of S along Z. Thus, we 
make the assumption that Z is compact. 


Theorem 1.12. The Wick expansion Wz{K) of I depends only on the diffeomorphism class 
of the integrand of I. 


Proof 


This is an automatic consequence of Lemma 1.9 □ 


Note that because of the above theorem, it was justifiable to speak of the Wick expansion 
as being a function of I instead of its integrand, since in case I is convergent, its numerical 
value depends only on the diffeomorphism class of its integrand as well. 

Altogether, we see the formal analogies between the Wick expansion and integration. 
Both are coordinate independent and satisfy analogous properties with regards to integra¬ 
tion by parts and fiber integration. We can summarize these observations in the following 
informal principle: 
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Formal Integration Principle: Any natural identity on integrals yield a corresponding 
identity on Wick expansions. 

The next section presents a further manifestation of this principle. 

2. The Wick Expansion and Gauge-Fixing 

We now consider integrands that are invariant under a (not necessarily compact) Lie 
group G of symmetries. Such a symmetry is regarded as a redundancy among the variables 
occuring in the integrand. In the context of integration, this means we should factor out 
the contribution of G to the associated integral. In the context of the Wick expansion, this 
means we need to eliminate the degeneracies in the action arising from G-invariance. For 
both these situations, the Faddeev-Popov procedure allows us to perform the required oper¬ 
ations just described in terms of a choice of gauge-fixing, which we define in a moment. This 
procedure actually has two different incarnations which we refer to as the slice and weighted 
Faddeev-Popov procedures. The distinction between these two methods and whether they 
should be applied in the context of integration or the Wick expansion is blurred in the 
literature, and our goal here is to give a clear and unified treatment that emphasizes their 
different features. This will pay great dividends later when we study the Faddeev-Popov 
procedure in the setting of quantum field theory. 

We first work in the context of integration. Moreover, we begin by describing the slice 
Faddeev-Popov method. The setup is as follows. We have a manifold M equipped with a 
volume form dV that is preserved under a left-action of G. We suppose G acts freely on M 
else we can restrict to a subgroup of G. Consider the integral 

1= [ dVfix) 

Jm 

where / is a G-invariant function. To eliminate the redundancy arising from G-invariance, 
we can choose a local slice S for the G-action, that is, a locally closed submanifold S C M 
which is transverse to the action of G (i.e. the G-orbit G • S through S is diffeomorphic to 
G X S). We refer to this procedure of choosing a local slice as gauge-fixing, since its analog 
in inhnite dimensions is what one does when choosing a gauge-fixing condition. 

Having chosen a local slice S, we can replace the part of the integral I over G • S with a 
suitably weighted integral over S. The most notable feature of this gauge-fixing procedure is 
the presence of a determinant. Such a determinant can be described implicitly using merely 
the fact that S is transverse to the G-action or else explicitly with the aid of a G-invariant 
Riemannian metric on M such that dV is the associated Riemannian volume form. This 
determinant takes into account that the volumes of G-orbits vary in M, so that a weight is 
needed when passing from an integral over M to that on S. (For the case of G noncompact, 
one instead considers the “ratio” of volumes of different G-orbits). Such a determinant is 
referred to as the Faddeev-Popov determinant. 

We first prove an “implicit” version of the slice Faddeev-Popov method and then derive 
an explicit Riemannian version. While the implicit version makes no auxiliary choices, it is 
the explicit version that is more useful in practice for doing computations. In what follows, 
we assume without loss of generality that M is globally a product 

M^GxS. 


( 2 . 1 ) 
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This is aiways true iocaiiy, and we can assume it is true giobaiiy in the context of integration 
by using a G-invariant partition of unity. We refer to S in (2.1) as being a giobai shce. 
Given a giobai siice S, the muitipiication map 


p : G X S ^ M 

{9,w) ^ 9-w, 


is a G-equivariant diffeomorphism, where G acts on G x 5 by ieft muitipiication on the 
first factor. Pick any ieft-invariant voiume form dVc on G. In doing so, we can define the 
voiume form on S via 

p*{dV) = dVc X dVg^^. 

Indeed, because p yieids an identification Tg-^M = TgG © T^S at every point g ■ w G M, 
an eiement in the top exterior power of T*,^M and T*G determines one for T^S. 

Theorem 2.1. (slice Faddeev-Popov formula, implicit version) For any two global slices S 
and S', we have 


[ dV^/^f{w) = [ dV^^/^fiw). 
Js Js' 

Furthermore, if G is compact, then 

[ dVf{x) = VoI(G) / dV^^^fiw). 
Jm Js 


( 2 . 2 ) 


(2.3) 


Proof The case G compact foiiows readiiy from the fact that we can integrate aiong 
the fibers of the projection tt : M —)• 5 and dV^^^ satisfies 


7r,(dP) = Voi(G)dP//^. 


Thus, (2.2) and (2.3) foiiows from ( |1.14 ). 

For generai G, we proceed as foiiows. Suppose we have two giobai siices S and S'. This 
means there exists a “gauge transformation” 9 : S ^ G that yieids for us an induced 
diffeomorphism 


0 : 5 ^ 5 ' 

w I—)• 9{w) ■ w 


reiating the siices via the G-action. We want to show that 

e*{dv^^^) = 

Define the map 


(2.4) 


(2.5) 


Le : M ^ M 

g ■ w g ■ w' = g6{w) ■ w. 
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We have the commutative diagram 


GxS 


M 


GxS' 


Lg 


M 


The vertical maps are volume-preserving by definition, since the volume forms on the prod¬ 
uct spaces are induced by pullback from multiplication. Thus, showing (2.5) is equivalent 
to showing that Lg is volume preserving. 

Given x G M, we have the map px ■ G ^ M, g ^ g ■ x determining the G-orbit through 
X. Define its differential to be 


I'x '■ Q ^ M 

ix{X) = ^ • X, X G 0. 

dt t=o 

giving via the infinitesimal action of G on M. In this way, each tangent space to the G-orbit 
through X can be identified with the space of right-invariant vector fields of G via the map 

L. 

We can write an arbitrary tangent vector u G Tg-^M in terms of its components with 
respect to the G-invariant distributions given by the tangent spaces to the G-orbit G • w 
and S, 

u = g*Lw{X) + g^(y), X G 0, u G r^5. 

Likewise, we can write a tangent vector at as {g6{w))^iu,{X) -|-((D^0)(u)). 

We have {Du,Q){v) = Rw{v) + were Rw maps v into the image of 6{w)^: o Lw 

Piecing all these decompositions together, the derivative of Lg{g ■ w) = g ■ Q{w) at g ■ w is 
given by 


D. 


^g-wLg ^g* iw (^)) 9 * 


({9d{w)g ^)^{g^Lw{-)) g*Rw{-) 

V 0 {gd{w)g-^)^{g^{-)) J \v J ' 


Being upper triangular, it determines the same map on top-degree forms as {g6{w)g~^)^:. 
The latter is volume-preserving since left multiplication by G is volume-preserving. Thus, 
Lg is volume preserving, thereby establishing (2.5). 

Hence, for / a G-invariant function, 



The disadvantage with the above formulation of the Faddeev-Popov procedure is that 
it does not express the volume form in a very explicit manner (especially for gen¬ 

eralization to quantum field theory). However, if we endow M with a G-invariant metric, 
with dV the associated G-invariant Riemannian volume, we can describe dV^^^ explicitly 
as follows. The Riemannian metric on M restricts to a Riemannian metric on S and so 
induces its own volume form dVs on S. 
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Definition 2.2. The Faddeev-Popov determinant Jsiw) is the function on S defined by 


dVs^i"^ = Js{w)dVs,n,, wgS. 


( 2 . 6 ) 


An explicit characterization of the Faddeev-Popov determinant, which is indicative of its 
name, is as follows. The G-invariant volume form dVc on G is determined by the volume 
form dVg it induces on its Lie algebra. For w G S, we can compose Lw with the orthogonal 
projection onto T^S, the orthogonal complement of T^S inside T^M. This yields for us 
the map : g —)• T^S. It is an isomorphism since S is transverse to the G-action. Let dVj- 
and dVs denote the volume forms on T^S and S, respectively, determined by the metric 
induced from M. Then we have 


{iindVlJ = Js{w)dV,, wgS. (2.7) 

Explicitly, let et be a basis for g* such that dVg = Ae* and let v*{x) be an orthonormal 
coframe for S, i.e., the v* {w) form an orthonormal basis of Ann(T^5) C T*M. In particular, 
d^s,w = Define Aij{w) by 


Then 


Js{w) = det Aij{w). 


( 2 . 8 ) 


An automatic consequence of Theorem 2.1 and the definition of Js{w) is the following 
Riemanian version: 


Theorem 2.3. (slice Faddeev-Popov formula, Riemannian version) For any two global 
slices S and S', we have 


J^dVsJs{vj)fiw) = 
Furthermore, if G is compact, then 


dVs'Js'{w)f{w) 


(2.9) 


dVf{x) = Yo\{G) [ dVsJs{w)f{w). 


( 2 . 10 ) 


JM JS 

Remark 2.4. In some formulations of the Faddeev-Popov formula, a square root of a 
determinant appears instead of a determinant. Namely, one has the formula 


Js{w) = detD2((i^)%^) ( 2 . 11 ) 

where (i^)* is the adjoint of : g —)■ T^S. The square root occurs because the volume 
form associated to a metric tensor gij receives a factor of (gij), and the pullback by 

of the metric on T^S is given by 


iyjU, iyjV 


T^S 



U,V G Q. 


Next, we discuss the weighted Faddeev-Popov procedure which is a generalization of the 
slice method. We discuss both these methods since there is a practical distinction between 
these two methods when applied to quantum field theory, as we discuss in Section 4.3 For 


the weighted method, instead of fixing the gauge by having an integral over S, we break 
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gauge-invariance by introducing a suitable weight function on M which is non-constant 
along gauge-orbits. The slice method is recovered by letting this weight function tend to 
the delta-current determined by the chosen slice. 

More precisely, let T : M —?■ Q be a gauge-fixing function from M to some target manifold 
Q such that 

(i) S = F“^((7o) with go a regular value (i.e. dFw '■ T^M Q is surjective for all 

w G 5); 

(ii) F restricted to some (and hence every) gauge-orbit has nonzero degree. 

Here, the degree deg(F), in the case of G and Q compact, is defined in the usual way: if 


F o pyj ■. G ^ Q 

g I—)• F{g ■ w) w G S 


( 2 . 12 ) 


pulls back a unit-volume top-degree form on Q to d times such a form on G, then deg(T’) = d 
is the degree. The independence of the choice of w follows from the fact that the degree is a 
homotopy invariant. One can also adapt the above setup to the situation of noncompact G 


by requiring (2.12) to be proper and working with compactly supported top-degree forms. 


but we will not do so here since we are ultimately interested in the Wick version for which 
G can be arbitrary. 

Pick a unit-volume top-degree form dq on Q and a weight function : Q —)• M such that 
f dq(p{q) = 1. The following tells us that by inserting into I the weight (p{F{x)) along with 
the determinant det(dF o defined via 

{dF o ix)*{dqF{x)) = det(dF o Lx)dVQ, 

we get a weighted integral independent of the choice of F and ip up to an overall constant: 


Theorem 2.5. (weighted Faddeev-Popov formula) We have 


IM 


dVifiFix)) detidF o L^)f{x) = deg(F) / dV^>^f{w) 


Proof Pull back the integral over M in (2.13) to 5 x M via p : G x S 


w G S, when we integrate over its G orbit, we have to perform the integral 

[ dVc ifiF o py,) det{dF o iyj) = f {F o py,)* {dq (f{q)) 

Jg Jg 


= deg(F) / dqip{q) 
JQ 

= deg(F). 


Peforming the remaining integral over S yields the result. □ 


(2.13) 
M. Fixing 


Remark 2.6. Let us take a closer look at how our formulation of the Faddeev-Popov 
formula relates to the usual informal expression involving delta functions, which is expressed 



as 


dVd{F{x)) det{dF o ix)f{x). 


(2.14) 
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By slight abuse of terminology, one also refers to det(d-Fo i^) as the Faddeev-Popov deter¬ 
minant, although the term 6{F(x)), which restricts the integral to S, contributes a factor 


coming from the differential of F. We can regard (2.14) as arising from (2.13) by letting 
ip become a delta function at qq. We can also compare (2.14) directly with (2.9) to see 


how to make a direct translation between these two formulations. Here, for concreteness, 
we suppose F : M ^ q, with S = F~^{0). In local coordinates, write the 1-form dF^ 
restricted to S in terms of an orthonormal conormal frame for S, i.e. dF^ = for 
some matrix Bj. Then 6{F(x)) acquires an inverse determinant factor (detHp“^, which 
cancels the corresponding factor of detH* arising in det{dF o t^). The overall determinant 


factor we obtain is (2.8). In this way, we see that transformation properties of the 6{F{x)) 


and det{dF o l^) conspire to make (2.14) independent of the choice F used to dehne the 
slice S. 

We now intertwine the Faddeev-Popov formula with the Wick expansion. Consider the 
integral 


/ = 


IM 


dVf{x)e-^^^^/^ 


(2.15) 


where S and / are both G-invariant. We want to obtain a Wick expansion of this integral 
about a critical subset of S. Such a Wick expansion requires a choice of gauge-hxing: we 
need to eliminate degenerate directions arising from the G-invariance of S{x) in order to 
obtain a splitting of S{x) into a nondegenerate quadratic part and a higher order interaction 
part. 

Note that every critical set Z of S{x) is G-invariant. We say that Z is Morse-Bott 
G-nondegenerate if on the quotient M/G, the set ZjG is Morse-Bott nondegenerate with 
respect to the induced function S{x). Given such a function S{x), we can dehne a Wick 


expansion of (2.15) by adapting the Faddeev-Popov procedure dehned for integration. We 


can adapt either the slice or the weighted procedures, but unlike in the case of integration, in 
which a degree term appears that distinguishes these methods, the Wick expansion version 
of these procedures yields the same output. 

Let us begin with the slice version. On a small G-invariant patch of M, instead of I we 
consider the integral 

Is = j dvf 

It is this integral that we can Wick expand using the methods of Section o Indeed, 

Zs ■■= Zr\S 

is a Morse-Bott nondegenerate level set for S restricted to S and we can choose a Morse-Bott 
bundle structure tt : —)■ Z 5 as before, where Zs is a tubular neighborhood of Zs inside 
S. Given any two gauge-hxing slices S and S' ^ we can always relate them via a unique map 
0:5 


S' determined by the G-action as in (2.4). The next theorem tells us that Wick 


hber integration along a gauge-hxed critical set Zs is covariant with respect to changes of 
gauge. 

Theorem 2.7. Let S and f be G-invariant functions and Z a Morse-Bott G-nondenerate 
submanifold. Pick a slice S and a fiber bundle structure vr for Zs- Let S' be any other slice 
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and let vr' be the fiber bundle structure for Zs' obtained from vr by the G-action relating S 
to S'. In other words, we have the commutative diagram 

Zs Zs' 


Zs — Zs' 

Then 

vr'r(V5') = 0*(vrf(V5)), (2.16) 


where V 5 and V 5 ' are the integrands of Is and Is', respectively. In particular, if Zs = Zs', 
then Tr'^iVs') = (V 5 ). 


Proof 


This is an automatic consequence of the definitions, 


Lemma 1.9, and (2.5). □ 


For Zs compact, we can apply Definition 1.11 to obtain the Wick expansion of Is about 
^5: 


Wzs{h) = 



(2.17) 


That is, the coefficients of tt^ (V 5 ) are top-degree forms on Zs, and we integrate along Zs 
to obtain a formal series in h. 


Theorem 2.8. The Wick expansion Wzg{h) is independent of the choice of gauge-fixing 
slice S. 


Proof 


Equation (2.16) shows that the integral (2.17) is independent of the choice of5. □ 


We may thus refer to Wz^ as the gauge-fixed Wick expansion of I about Z with respect 
to 5. It is a well-defined formal series in h, which requires a choice of gauge-fixing condition 
S for its construction but is independent of the choice made. It is built out of Wick 
contractions applied to the gauge-fixed integrand Vs, and moreover its gauge-invariance 
arises from the covariance properties of Wick fiber integration. This involves formal algebra 
along the fibers rather than integration. Hence, the gauge-fixed Wick expansion serves as 
the correct finite-dimensional analog of gauge-fixed perturbative path integrals, as we will 
discuss in Section 4. 

Having explained the Wick version of the slice Faddeev-Popov procedure, we now turn 
to the Wick version of the weighted Faddeev-Popov prcoedure, whose ansatz turns out to 
be more useful for the perturbative quantization of gauge theories. Here, in augmenting the 
integrand of I, we pick a gauge-fixing function T : M —)• Q and weight function y? : Q —)• M 
as before, with S = F“^(go). We also make the additional assumption that (p is of the form 
where h{x) has a nondegenerate critical point at go and Cf, is defined to be the 
Wick expansion of f dqe~^^‘^'t/^ (thus is a Laurent series in h}/'^). 

We thus consider the weighted integral 

4 ,f = J dVip{F{x))det{dF o i,,)f{x)e-^^^^/\ 
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It has Morse-Bott nondegenerate set Zs since ip{F{x)) modifies the action S{x) to a new 
one that is nondegenerate in the G-directions. 


Theorem 2.9. The Wick expansions of and Is about Zs agree. In fact, if n is a fiber 
bundle structure for Zs inside of M (as in Definitional^, then 

(^dV(p{F{x))diei{dFoix)f{x)e~^^^^^^'^ (2-18) 

is independent of the choice of F and ip satisfying the requisite hypotheses. Consequently, 
the Wick expansion of Ip^^ is also independent of F and p. 


Proof We know that the Wick expansion is independent of the choice of coordinates. 
So pull back the integrand of I^p^p by p : G x 5 —)■ M. Since h{F{x)) is nondegenerate 
along the G-directions and constant along S, the Wick expansion in the G x S coordinate 
system factorizes in the G and S directions. When we perform a Wick expansion in the 
G-direction, we have to Wick expand 


c 


-1 

h 



dVc det(dF o Ug.^) 


about w € S. By coordinate-independence {g ^ F{g ■ w) is a diffeomorphism in a neigh¬ 
borhood of the identity element of G) this equals the Wick expansion of 


cl^ [ dqe-^^^'^!^ 

Jq 

which is 1. The leftover Wick expansion involves Is. D 


Thus, the Wick version of the weighted Faddeev-Popov procedure only sees a “local” 
degree (equal to one) instead of the global degree that occurs in the integration version. 

3. The Wick Expansion and Integral Asymptotics 

In the previous sections, we considered the Wick expansion of an integral as a purely 
formal series in h. This allowed us to consider properties of Wick expansions independently 
of the convergence properties of such integrals. In this section, we consider those integrals 
which are convergent for h along a ray in the complex plane. This makes 


I{h) = / d^xf{x)e-^^^^/^\ 


(3.1) 


a function of h. We consider h along different rays because the cases h real and imaginary 
correspond to integrals arising from Euclidean and Lorentzian physics, respectively. Inter¬ 
mediate cases correspond to integrals obtained through Wick rotation. In what follows, we 
always assume S{x) is real-valued. 

As is well known, the Wick expansion provides an asymptotic expansion for I{h) as 
h —)• 0, in the sense of (3.3), under suitable hypotheses. In this context, the method by 


which it is shown that the Wick expansion provides the correct asymptotics goes by several 
names, including saddle point approximation, stationary phase, or steepest descent. The 
asymptotic nature of the Wick expansion makes it unsurprising that the Wick expansion 
has the covariance properties that it has. Indeed, the Wick expansion automatically inherits 
these properties from the corresponding ones for ordinary integrals when the former is an 
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asymptotic expansion of the latter. However, what the previous sections show is that, 


essentially as a result of Lemma 1.4, the covariance properties of the Wick expansion, with 
h a formal parameter, are of a purely algebraic nature and so hold without any additional 
hypotheses on the underlying integrand. 

For completeness, we provide a proof of the asymptotic properties of the Wick expansion 
in the lemma below. We treat the case of nondegenerate critical points, everything being 
exactly analogous for the Morse-Bott case. First, we deal with the case when f{x) has 
compact support and we work locally near a critical point xq of S{x). This is the most 
natural case, since for oscillatory integrals in which h is imaginary, the term e-‘S'(x)/;i 
not decay at infinity and thus one has to truncate the integral smoothly. 

Write the Wick expansion as 




Lemma 3.1. Consider the integral 


fc=o 


I{h) = / 


(3.2) 


with f{x) compactly supported. Suppose S{x) has a unique nondegenerate critical point xq 
on the support of f{x). Let /i —)• 0 along a ray in the complex plane, where ifReh > 0, 
suppose further that xq is a minimum of S{x). Then Wxoih) is an asymptotic series for 
I{h) in the sense that for every N > 0, 


I{h) - Yj 


(3.3) 


k<N 


Proof Without loss of generality, we can suppose xq = 0. The asymptotics of I{h) are 
unchanged if we replace I (h) with 

I{h) = J d‘^xif{x)f{x)e-^^^'>/^ ( 3 . 4 ) 

where f’ix) is any bump function that is identically one in some neighborhood of the origin. 
For Re/i > 0, this is because for x in the support of /(x) and bounded away from xq, 
Q-S{x)/h ^ ^-S{xQ)lh^-Clh £ qj . gQjj^g positive constant C. For Reh = 0, we use the standard 
integration by parts trick. Starting with d = 1, we have 

j dx[l-fj{x)]f{x)e~^^^'>/^ = j dx[l-'ip{x)]f{x)(^-^l^-^^ ( 3 . 5 ) 

for arbitrary N . On the support of [1 — ?/;(x)]/(x), we have that ^'(x) is nonzero, in which 
case integration by parts shows that the above integral is 0{h^). Since N was arbitrary, 
such a term is asymptotically zero. For d > 2, one can adapt this technique to draw the 
same conclusion, see O Chapter VIII]. 

The inverse function theorem allows us to choose local coordinates near the origin in which 
5(x) = ^H(x,x). By making the support of ip{x) as small as we like, we can suppose that 
such coordinates have been choosen globally. This is a more convenient coordinate system 
to establish that the asymptotics of the function I{h) is equal to the Wick expansion of I{h) 
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regarded as a formal integral. Since the Wick expansion is independent of the coordinate 
system chosen, our choice of coordinates does not affect the result. 

Pick (5 > 0. Write 

e^^^^"f{x) = PN{x) + RN{x) 

where Pn is the degree N Taylor polynomial of centered at the origin and Rn 

the remainder. So then 

J d^x^{x)f{x)e-^^^^/^ = J 

+ J d^xi^{x)RN{x)e-^^^^\-^^^’^'>/‘^^. (3.6) 


For the first term, we apply Wick’s theorem using the quadratic form 2(5|xp + A{x,x)/h, 
which has positive definite real part for 6 small. Sending (5 —)• 0, we obtain leading terms 
of the Wick expansion. For the second term, if Reh > 0, we have pointwise exponential 
decay of the integrand. If Reh = 0, the same estimate used t o co ntrol (3.5) shows that 
the integral is 0{hM) for all M > 0. Finally, the third term of (3.6) is For 

Reh > 0, this follows from the rescaling x i—)• (Reh)^/^x and using Rn{x) = x^^^r]{x) with 
r] smooth. For Re h = 0, one has to work harder, see m p. 335]. Since N was arbitrary, 
this establishes the lemma. □ 


Next, we want to remove the hypothesis that f{x) is compactly supported. Additionally, 
we want to impose conditions such that the sum of Wick expansion of I{h) at all of its 
critical points yields the full asymptotics of I{h). For simplicity, we assume S{x) has 
finitely many critical points Xc and that all of them are nondegenerate. (The asymptotics 
of I{h) about degenerate critical points can be studied, but they cannot be analyzed via the 
Wick expansion.) Moreover, we assume S{x) —>• oo sufficiently rapidly at inhnity. Thus, 
for Reh > 0, as long as /(x) does not grow too quickly, (3.1) is well-defined. For Reh = 0, 
the function is no longer damping at inhnity but highly oscillatory. Nevertheless, 

so long as S{x) is suitably well-behaved (depending on f{x)), we can make sense of the 
regulated integral 


1(h) = lim Ie(h) 

e-5-O 

=: lim [ d'^xijj{ex)f{x)e~^^^'^/^. 

e ->0 / 


(3.7) 


Namely, we need ^ 7 ^ and all its derivatives to decay sufficiently rapidly at inhnity relative 
to /(x), so that we can use integration by parts to control (3.7) uniformly in e. 


Theorem 3.2. Assume the above hypotheses and those of Lemma 3.1, Let {xc} denote 


the set of critical points of S{x) z/Reh = 0 or else the set of minima of S{x) z/Reh > 0. 
Assume all the Xc are nondegenerate. Then Wx^{h) is an asymptotic series for L(h). 
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Proof The assumed decay properties allow us to control the final two terms of (|3.6) 
as before, where we replace the function with '0(ex) and let e 
0 (ji{d+N+i)/ 2 ) uniformly as e —>■ 0. □ 


0. These terms are 


4. Remarks on Quantum Field Theory 


We now apply our previous mathematical analysis to the setting of quantum field theory. 
We assume the reader has experience with quantum field theory so that we can be allowed 
to provide succinct commentary rather than a self-contained exposition. First, we review 
how the infinite-dimensional Wick expansion works in the setting of quantum field theory. 
After this, we carry out our main goal of explaining how a proper distinction between the 
Wick expansion and integration (as inspired by our previous finite-dimensional analysis) 
provides a conceptually sound treatment of formal path integral manipulations. 

In quantum field theory, one wishes to apply the Wick expansion to a path integral, i.e., 
an integral over the space of field configuration^ The path integral is formed out of a 
classical action 5, from which we obtain the classical equations of motion by considering 
those configurations for which the action is stationary. One may then attempt to Wick 
expand the path integral about these stationary configurations by a formal application of 


Definition l.c * Aside from an overall normalization, this is possible as long as we can 
make sense of the terms appearing in the Wick formula (1.5). The problem of course is 
that now the right-hand side of (1.5) is an integral of a product of distributions, since the 
inverse matrix is to be replaced with an appropriate Green’s function (the propagator) 
and the sum over indices become operator multiplication and integration. These are the 
familiar divergences appearing in Feynman diagrams. Thus, making sense of the Wick 
expansion, and indeed all other quantities in quantum field theory, require the proper use 
of regularization and renormalization. 

There are three common methods one can use to regulate path integrals: 

• lattice regularization: replace the continuum theory with one defined on a (finite) 
lattice, so that the path integral becomes a product of ordinary integrals; 

• propagator regularization: replace the integral kernel of the propagator with a 
smoothed out version, using e.g. a momentum cutoff or a heat kernel regulator; 

• dimensional regularization: do not modify the propagator, but analytically continue 
the linear operator f d'^x for integer d to complex d (a rigorous treatment can be 
found in [sQ. 

In all of these regularization schemes, there is a regulatory parameter e that serves as an 
ultraviolet regulator which is removed as e —0. (One may also need an infrared regulator. 


^Properly speaking, this terminology is appropriate for Lorenztian theories in which there is a notion of 
time, hence the notion of a path. We use this terminology to include Euclidean theories. 

o 

One almost always has a moduli of solutions to the classical equation of motion, whose finite-dimensional 
analog yields the Morse-Bott situation of Section fl.l[ We apply Definition |1.3| to compute Feynman diagrams 
evaluated on on-shell configurations, i.e. those satisfying the equations of motion, which is analogous to 
performing Wick hber integration (Feynman diagrams) and then considering Z (on-shell configurations). 

“^Note that in , the dimensionally regularized integration operator J ddx is only a linear operator when 
restricted to a suitable space of functions, e.g. rational functions. This is sufficient for the kinds of integrals 
that appear in translation-invariant field theories, whose integrands have Fourier transform equal to rational 
functions of momentum variables. 
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which can always be fixed as one removes the ultraviolet cutoff.) For instance, e can serve as 
a lattice spacing, the inverse of the momentum cutoff, or else one might work in dimension 
d — e. With a regularization scheme in place, this renders all integrals occurring in the 
regulated Wick expansion Wf:{h) finite for e > 0. Thus, We(h) becomes a well-defined 


formal series in h since the Wick formula (1.5) generates well-defined coefficients at every 
order in h, i.e. every loop order (modulo an overall normalization constant arising from a 
determinant, which we ignore from now on). 

Next, renormalization begins by making the action depend on e, in which we write = 
S + CT^, a sum of the original action plus counterterms CT^. Such counterterms, which 
diverge as e —)• 0, are arranged so that the renormalized Wick expansion (the Wick expansion 
associated to S^) has an e —)• 0 limit as a formal power series in h. 

This is the general schema by which the standard path integral approach to perturbative 
quantum field theory proceeds. Each theory requires the appropriate technical execution 
of the above procedures, along with any additional requirements imposed by other consid¬ 
erations such as symmetry. The successful execution of such steps provides a perturbative 
definition of a path integral, one which a priori is disconnected from trying to perform inte¬ 
gration. From this, the basic principle behind making sense of path integral manipulations, 
in direct analogy of the Formal Integration Principle of Section 1, can be summarized as 
follows: 

Formal Path Integration Principle: Any formal manipulation of a perturbatively de¬ 
fined path integral yields provisional identities whose legitimaey depend upon an analysis of 
the regularization and renormalization scheme employed. 


While those knowledgeable in quantum field theory are well-aware of this fact, it is unfor¬ 
tunate that this tacit philosophy is often obscured by formal notation or else interwoven 
with other procedures that have no rigorous analog. The approach we adopt in this paper 
yields a most natural elucidation of the Formal Path Integration Principle, since it cleanly 
separates the distinction between integration and formal integration already present in fi¬ 
nite dimensions from the technicalities of regularization and renormalization occurring in 
infinite dimensions. 

We now investigate how a variety of path integral manipulations can be interpreted 
rigorously by conscientiously using the Wick expansion together with the above Formal Path 
Integration Principle. We proceed as follows. Section 4.1 provides a warmup that shows 


how the different regularization schemes we outlined above work in generating the Wick 
expansion. Section [4.2| shows how the subtleties involved in understanding Ward identities 
can already be gleaned from the finite-dimensional distinction between integration and the 
Wick expansion. In Sections 4.3 and 4.4 the theorems we proved in Section are brought 
to their full light in terms of interpreting perturbative path integrals through the Wick 
expansion and not integration, and we explicitly mention some of the shortcomings in the 
literature. 


4.1. Change of variables in the path integral. For simplicity, we illustrate how to 
make sense of changes of variables in the path integral for the case of scalar field theories on 
M'^, though everything we discuss readily carries over to more complicated theories. That 
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is, we consider path integrals of the form 


where 


1= ^ f 


S{(j)) = I d^’-x ^(p{x)A4>{x)+ + V{4>{x)) 


(4.1) 


V and O are some polydifferential functions of 4>, and Z is a normalization constant. What 
follows is a toy computation to make explicit how regularization and the Wick expansion 
of path integrals intertwine. We then make comments about the more general situation 
afterwards. 

To begin, consider a free massless theory with action 


S(4>) = I 


d'^x 4>{x)A(f>{x). 


With 0 = 1, the Wick expansion of (4.1) is trivial (i.e. is a constant which we normalize to 
one) since there are no interaction terms. However, suppose we make a change of variables 
4>{x) I—)■ 4>{x) + (p^{x). The new action will no longer be a free theory and so will generate 
Feynman diagrams via a Wick expansion of the path integral about the zero field. 

We can gain insight by seeing what happens in the one-dimensional case. The substitution 
X ^ X + x^ does the following: 

1 f , _.i.2 /ofc 1 


dxe-^ 


dx(l + 


V^nh J \/27rh 

By Theorem |1.5[ the Wick expansions of the two integrals about xq = 0 are identically 
the same (and in this case, the integrals are themselves convergent). But while the Wick 
expansion of the left-hand side is trivial, the right hand-side involves cancellations of di¬ 
agrams to all orders in h, since we have interaction terms. For instance, at first order 
in h, the Wick expansion receives a contribution from and This yields 3h and 

—3h, respectively, which cancel. Likewise, at order h?, we have to apply Wick’s theorem to 

-3x2 .1^ + 1 - iJ, which yields ff{-3 • 5!! + ^7!! - ^5!!) = 0. 

In quantum field theory, the combinatorial factors occurring in the Wick expansion are the 


same as those that appear in the Wick formula (1.5). It is the Wick contractions appearing 


in (1.5) that are affected by which regulator we choose. We illustrate what happens using 


all three methods previously discussed. 

For lattice regularization, we have lattice points Xi belonging to a finite-size box embedded 
within a rectilinear lattice with lattice spacing 6x. Let A denote the lattice Laplacian with 
periodic boundary conditions (so that A = Ajj is a matrix) and let (•, •) denote the lattice 
inner product defined by 

if,9) = '^{^xYf{xi)g{xi). 


The lattice partition function is then 

Z= j 











PATH INTEGRALS 


25 


where d(p{xi) is Lebesgue measure for the real variable (p{xi). Under the change of coordi¬ 
nates (l){xi) I—)• the integral changes to 

i 

Let Gij be the inverse matrix of Ajj restricted to the orthogonal complement of constant 
functions. In the continuum limit dx —)■ 0, we have 

^(5x)‘^/(xi) ^ j (fxf{x), 

so that Gij{5x)~^ is the lattice version of the continuum Green’s function: 

'^Gijf{xj) = '^{5xfGij{Sx)~‘^f{xj)^ j (fxG{x\x)f{x). 
j j 

Moreover, Gij{6x)~'^ is the matrix we use to perform Wick contractions when we Wick 
expand. 

At order h, we have to consider Wick contractions of 3(p‘^{xi) and the quartic interac¬ 
tion. Wick contraction of yields 3Y^^Gii{5x)~'^. It is a regulated version of its 

divergent counterpart in the continuum limit, which is given by 

3((5x)-‘^ V((5x)‘^Gii((Ix)-'’'^ 3(5‘^(0) [ d^xG{x,x) (4.2) 

. Jm 

where the latter integral is performed over the torus M given by the periodic identification 
of the box defining our lattice. Wick contraction of the quartic interaction yields 

- ^((5x)‘’*(/>^(xi)A(/>(xi) —^ -3'^{6x)'^Gii{6x)~‘^AijGji{6x)~‘^ (4.3) 

i i,j 

= -3Y,Gii{Sx)-‘^. (4.4) 

i 


This exactly cancels out the previous diagram. 

For propagator and dimensional regularization, one starts with the continuum theory. 
The path integral transforms as 

J D(l) ^l + ,5('^)(0) J d^x3(j)'^{x)^ g-[i(<^,A<A)+(0,A<A3)+|(0LA03)]/?i^ 

One can justify this formal manipulation on the basis of the formal continuum limit of the 
lattice theory above. The difference is that previously, all quantities were well-defined on a 
finite lattice whereas now one has to “undo” these formal manipulations with a regulariza¬ 
tion scheme different from the lattice. 

For propagator regularization, this means regulating 


G{x,y) = 


1 

(2^ 

1 

( 2 ^ 


eha^-2/) J_ 

|p|2 

d'^p 


(4.5) 

(4.6) 
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by taming the integrands in these expressions. In the present situation, the two Feynman 
integrals we have to consider are 35^'^\0) f d^x G{x, x) and —3 f d^x G{x, x)AG{x, x) arising 
from Wick contracting f d‘^x34>{x)‘^ and — f d^x /2h, respectively. Formally, 

5('i)(0)-AG(x,x) = 0 (4.7) 


and so the two Feynman integrals cancel. To make (4.7) more meaningful, we can regulate 


the theory in such a way that if we replace (5^'^^(0) and G in (4.7) with regulated versions. 


call them and G^, we still get exact cancellation. More generally, what we need is that 

lirn ^(5^‘^^(0) — AGe(x, = 0. (4.8) 

Such a regularization procedure makes rigorous the cancellation of two infinite quantities. 


integrands of (4.5) and (4.6). 


It is straightforward to obtain (4.8) by inserting the appropriate momentum cutoff in the 


In dimensional regularization, one obtains a Wick expansion using the integration oper¬ 
ator f d'^x with d complex. In this setting, it turns out there are no power divergences, 
i.e. f d^p\p\‘^°‘ = 0 for all real a. Thus, both and AG{x,x) are zero individually, 

i.e. both diagrams under consideration vanish. Strange as this may seem, so long as f d'^x 
is a consistenlj^ linear operation, this provides one with a well-defined rule for performing 
the Wick expansion. While the author has personal reservations about using dimensional 
regularization due to this strange elimination of divergences, it is the presence of such for- 
tuitious eliminations of divergences which makes dimensional regularization populaij^ 


Having discussed the above example, it is no more difficult conceptually to consider 
transformations of the form (f> (p + F{(p) with F a polynomial that is at least quadratic. 
Such a transformation is invertible as a power series in </>. One can inspect the diagrams 
one generates when performing such a change of variables and inspect how regularization 
meshes with the terms of the Wick expansion. Such an inspection has to be done on a case 
by case basis, with formal manipulations of the path integral being a way of delaying, or in 
many cases, masking this inspection. One could also consider transformations of </> in which 
the linear term is not the identity, but some general linear transformation L(p. In the Wick 
expansion, one picks up a determinant factor from this linear transformation. This must 
also be regulated appropriately. 

4.2. Ward identities and Schwinger-Dyson equations. One often wants to show that 
the path integral of a Lie derivative is zero, where the Lie derivative arises from some 
infinitesimal symmetry or infinitesimal change of variables. The resulting identities are 
referred to by various names depending on their context, which we recall in Table 1. Despite 
the many names, it is customary to refer to the general collection of identities obtained 
through infinitesimal variations of the path integral as simply Ward identities. 

®See also for some difficulties with dimensional regularization if not done properly. 

®In addition, dimensional regularization preserves gauge symmetry, which makes it a standard choice 
for quantizing gauge theories. We should note that dimensional regularization, at least as presented here, 
only applies to translation-invariant theories and thus has rather limited applicability. In contrast, lattice 
regularization and propagator regulation via heat kernel methods are robust and can be applied to general 
situations. For some comparisons between different regularization schemes in specific models, see e.g. HIS]. 
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Identity 

Operation 

Schwinger-Dyson equations 

Ward identity 

Slavnov-Taylor identities 

varying a scalar field by translation (/)(x) i—)• 4>{x) + e{x) 

varying the electron by its phase T —?■ in quantum 

electrodynamics 

varying all fields by the BRST operator in gauge-theories 


Table 1. Various identities occurring in quantum field theory. 


For the case of Schwinger-Dyson equations, the identity we obtain is usually expressed 
in the form 


{T'' 


l = j{ o{^) 


SS{<P) \ 
H /’ 


(4.9) 


for O a general observable, where (O) is defined by (4.1). This equation is interpreted as a 
“quantum equation of motion”, since f or ex ample taking 0=1, this equation says that the 
expectation of is zero. Equation (4.9) is obtained by declaring (i) the path integral is 
invariant under the translation change of variables 4>{x) i—)• 4>{x) + e(x); (ii) the “measure” 
D(j) is invariant under translation. Differentiating with respect to e at e = 0 yields 

(4.10) 

(4.11) 

(4.12) 


= ^ 6 


0 = 61 
1 
Z 
1 

Z 


.-sW/h 


D(j) 


(scm 
V H 


— h ^0{(p )—I e 


.-sw/n 

H ) 


We see that (4.9) is a formal analogue of (1.8). As such, (4.9) has to be interpreted in the 
sense of formal power series in h, in which each side generates its own Wick expansion. Thus, 
equation (4.9), which holds formally, holds legitimately if one can arrange for regularization 
and renormalization (as done in the previous section) in the Wick expansion to be carried 
out so as to render (4.9) true. For instance, using lattice regularization, (i) and (ii) are 
obviously true since Lebesgue measure at each lattice site is translation invariant. In fact, 
there will in general be correction terms to (4.9) due to the fact that nonlinear functionals 
of the field need to be renormalized. 

It is important to note that many treatments drop the explicit dependence on h in (4.9), 
treating h as a number (such as i [161 Ch 9.6]). This blurs the formal power series nature 
of the involved quantities. Moreover, we also have to remember that the Wick expansion 
is usually about a moduli of configurations. In the finite-dimensional situation of Section 
2, the analogue of Lemma j 1. 4| would hold only if one could determine that the integral of a 
total derivative on the on-shell space Z vanishes as well. This is also blurred in the formal 
notation (4.9). Any rigorous interpretation of the provisional equation (4.9) must take these 
considerations into accounlQ 


•7 

There are algebraic approaches to perturbative quantum held theory with a cohomological emphasis, 
in which one considers algebras of observables equipped with a differential (see e.g. O HO]). With the 
appropriate setup in place, one regards a quantum expectation as taking the cohomology of this differential. 
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The other Ward identities can be handled in a conceptually similar manneij^ The fact 
that one varies by a symmetry of the action means that one expects to obtain identities 
without terms proportional to the equation of motion. However, if the “measure” Dx on the 
space of fields x is not invariant with respect to the symmetry, one obtains an anomaly which 
is a correction term appearing as a Lie derivative of the measure. In the finite-dimensional 
Wick expansion, such a Lie derivative can be computed using the algebraic rule Cy = diy, 
in which one obtains a divergence factor diV^{x) associated toV = V^{x)dl.. Indeed, 

Cx{dx^ ■ ■ ■ = diV'‘{x)dx^ ■ ■ ■ dx'^. 


It is the appearance of this divergence factor (suitably regulated) that can be generalized to 
the infinite-dimensional setting without requiring the existence of a measure for integration. 

The well-known chiral anomaly of QED falls under this analysis, in which an anomaly 
term arises from the local index density of a chiral Dirac operator |16[ Ch 19.2]. A general 
framework of how to properly treat the change in the “measure” Dx under a symmetry in 
continuum perturbative quantum field theory is provided by the Batalin-Vilkovisky formal¬ 
ism as developed in [6] . For a concrete example of this formalism in the context of nonlinear 
sigma models, see [T5] . 

Altogether, we see how various subtleties involved with obtaining and interpreting Ward 
identities can be inferred from the algebraic features of the Wick expansion. Trying to 
discuss these issues in terms of integration causes difficulties when trying to justify the 
integration of an total derivative being zero in (4.11) or the existence of a measure that 
yields an anomaly. 


4.3. Gauge-fixing of Feynman amplitudes and path integrals. As an application 
of our analysis of the Wick expansion in Section 2, let us provide some insights into the 
gauge-fixing procedures done in quantum field theory. 

Our first goal is to provide a hnite-dimensional interpretation of the gauge-invariance of 
on-shell Feynman amplitudes. Accomplishing this was in fact, one of our motivations for 
this paper, since conventional path integral methods obscure why the on-shell condition 
is crucial. The specific example we have in mind arises from quantum electrodynamics 
(QED). Here, there is a natural family of gauge-fixing conditions parametrized by ^ > 0 
with corresponding propagator 





(1 -g)pMp^ \ 

p2 J ■ 


(4.13) 


In the finite-dimensional setting, this would be analogous to replacing integration over the on-shell space Z 
in Definition 1 1.11 1 with passing to cohomology. This automatically makes integration by parts valid, i.e., the 
expectation of an exact term is zero, without any hypotheses on Z. However, such a cohomological setup 
has the disadvantage of losing connection with integration, even in the rigorous finite-dimensional situation. 

®For those theories involving ghosts and odd symmetries, to accurately mimic those theories in the finite¬ 
dimensional setting, one would have to invoke a version of the Wick expansion which includes fermionic 
variables. However, as fermionic integration is already purely algebraic, one can readily extend our finite¬ 
dimensional analysis of the bosonic Wick expansion to the case of supermanifolds by adapting the differento- 
geometric framework of supermanifolds in [S]. 
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Inspection shows that it is the Wick version of the weighted Faddeev-Popov procedure 
(adapted to quantum field theory) that is used to obtain the propagator (4.13). Indeed, in- 


1 

2^h 




serting the weight exp ( — f 
in ( |4l3 ). _ 

gauge-fixing functional is F{A^) = and varying ^ corresponds to varying the weight 


into the QED path integral introduces the ^-term 
Mimicking Theorem 2.9 in this infinite-dimensional context, the corresponding 


</?(•) = e I'l Theorem 2.9 shows that Wick fiber integration (generating Feynman di¬ 


agrams) and then restriction to the gauge-fixed critical configurations Zs of the action S 
(on-shell evaluation)|^ is independent of the variation of (/?. Of course, as with the previous 
section on changes of variables, one always has to inspect how regularization and renor¬ 
malization affect the translation of any result about finite-dimensional Wick expansions to 
infinite dimensional ones. Since our goal here is to highlight the conceptual fitness of the 
Wick expansion, we leave the details to the industrious reader. 

By contrast, the conventional description of gauge-fixing which proceeds by way of anal¬ 
ogy with the finite-dimensional integration version of the Faddeev-Popov procedure is mis¬ 
leading. Indeed, the only finite-dimensional gauge-invariant regularization of gauge theories 
(to the author’s knowledge) involves working with group-valued instead of Lie-algebra val¬ 
ued fields on a lattice. Since in most instances, the continuum limit of the lattice formulation 
has not been rigorously shown to approach the continuum theory, it is unclear how the in¬ 
tegration Faddeev-Popov procedure on the lattice is related to the Wick Faddeev-Popov 
procedure in the continuum [18]. Next, the Feynman diagrams generated from a gauge- 
fixed action do not depend on topological features of the group of gauge transformations 
or of the gauge-fixing function F ; the diagrams are generated purely algebraically from the 
gauge-fixed action. In finite dimensions, we saw that such a feature is true of the Wick 
version of the weighted Faddeev-Popov procedure but not of the integration version for 
which a topological degree term appears. Hence, once more it is the Wick expansion, not 
integration, that is most conceptually fit for describing how gauge-fixed path integrals are 
perturbatively evaluated. 

On the other hand, there are contexts in which the integration version of the Faddeev- 
Popov procedure is relevant to being able to evaluate path integrals nonperturbatively. 
Here, we have in mind the computation of the path integral formed out of the Polyakov 
action in string theory (as performed in [l2|). Whereas the Wick version of the weighted 
procedure is used to yield well-defined Feynman rules in perturbation theory as above, the 
integration version of the slice procedure (Theorem 2.3) rewrites a path integral on a large 
space in terms of an integral on a much smaller space using quantities that one can regu¬ 
larize, namely the Faddeev-Popov determinant as defined by (|2.8|) or (2.11). In particular. 


if M denotes the space of Riemannian metrics on a surface of genus h, and G is the 
group of diffeomorphisms and conformal rescalings, then M/G is the (finite-dimensional) 
moduli space of genus h Riemann surfaces, and the slice Faddeev-Popov method converts 
the ill-defined path integral over M to a sensible (nonperturbatively defined) integral over 
a slice S. Here, we can see why the Riemannian version of the slice method (Theorem 2.3) 


®Wick fiber integration yields a top-degree form on Zs, but we can identify it with a function on Zs 
once we choose a fixed reference G-invariant volume on Z. The Taylor series of the resulting function yields 
polynomials whose evaluation on Zs correspond to evaluating Feynman diagrams with elements of Zs placed 
on external legs. 
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is more useful than the implicit version (Theorem |2.1[ ). The latter directly involves volume 
forms (e.g. Lebesgue measure), which do not always generalize to the infinite-dimensional 
setting. But for the former, since infinite-dimensional spaces can still be endowed with met¬ 
rics, the corresponding Faddeev-Popov determinant can be assigned meaning via a suitable 
regularization, which allows for an evaluation of the resulting gauge-fixed path integral. 

Altogether, we see how the distinct Wick and integration versions of the Faddeev-Popov 
procedure in finite-dimensions provide differing insights into how to proceed in quantum 
field theoretic contexts. 


4.4. Eliminating fields by their equation of motion. Often times, a field only appears 
quadratically in a path integral and one wishes to “integrate out” this field. Two well-known 
instances of this are the Nakanishi-Lautrup field h occurring in BRST gauge-fixing and the 
B-field in BE-theories. In the first instance, one wishes to perform the path integral 




(4.14) 


so as to insert the gauge-fixing term in the path integral over A“ (we ignore 

overall constants in front of path integrals in what follows). In the second instance, one 
wishes to convert Yang-Mills theory to a first order theory by writing 


dAe 2 n f = J dAdBe^'^^ 


fll'' ' flU' J ^ flu' ' ^ IIU J 


(4.15) 


The usual physics justification for these procedures is that one is either performing Gaussian 
integration or else eliminating a field using its equation of motion: 


h- ^ -d^Al/ih 

The latter justification is a notational shortcut for the former justification, which in turn 
is improper from the point of view of honest integration, since the integrals over /i“ and 
involve the exponential of a positive definite (^ > 0) or else imaginary (in Lorentzian 
signature) quantity. Hence, the correct justification, comes from acknowledging that these 
procedures only involve the Wick expansion, which is well-defined regardless of the signature 
of the quadratic part of the action. 

We can abstract the above analysis to the finite-dimensional case as follows. Let T : A —)• 
Y be an invertible linear map between inner product spaces of dimension d, and V : X ^ Y 
a polynomial map that is at least quadratic. We have the action 

A(a:) = ;^|Ta: + Y(a:)|2 (4.16) 

2c 

where | • | is the norm associated to the inner product (•, •) on Y. One can Wick expand 
f about xq = 0. On the other hand, we can introduce an auxiliary variable 

y ^ Y and consider an equivalent “first-order formulation” (we think of T as being a hrst 
order operator). 

Namely, consider the first order action 

Spoix.y) = - (I + {Fx + V{x),y)^ 


(4.17) 
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The quadratic part of the action — ^ (c (?/, y) + 2 (Tx, y) ) can be identified with the matrix 


A = - 


0 T*\ 

\T c) 

via the inner product on X and Y. Its inverse is given by 




0 


1-1 


c{T*T)-^ 

This yields a complicated set of Feynman rules if we Wick expand Spo about (x, y) = 0. On 
the other hand, we know that the Wick expansion is invariant under changes of coordinates. 
Since 

- (I {y,y) + {Tx,y)^ = |y+ rx/cp + 

the change of variables y ^ y — Txjc diagonalizes the quadratic part of S{x,y). Further¬ 
more, being upper triangular, it is also a volume-preserving transformation. Thus the Wick 
expansion of SFo{x,y) is equivalent to the Wick expansion of 

S{x,y) = -||y|2 - {V{x),y) + ^(|Tx|2 + 2 {T{x),V{x))) 

Since the quadratic part of S{x, y) is diagonal, the Wick contractions of the x and y variables 
decouple. Wick expanding with respect to the y variables first, we obtain 


271/1 




(4.19) 


Thus, we recover the Wick expansion of S{x) up to an overall normalization. (In the first 
line, we regard {V (x), y) as an interaction. While this is not at least cubic in y, we can apply 
Wick’s theorem and get a well-defined series expansion in the formal variable h and in the 
x-variable, since V{x) is at least quadratic and so has no constant term.) Notice however, 
that (4.18)-(4.19) can only be interpreted as a formal Wick expansion. If Rec < 0, then 


(4.19) has the wrong sign for doing a Gaussian integral; if Rec > 0, then the integration 


over y in (4.18) has the wrong sign. In other words, the sign of c is always incompatible 


with either the y or x integral being convergent. It is unfortunate however that standard 
treatments of “integrating out” auxiliary fields describe this step (or some variation of it) 
as an honest integration procedure rather than as a Wick expansion. 

For a concrete instance, we can examine equation (15.7.5) from Weinberg’s textbook |20j . 
in which the gauge-fixing term exp(—^ f is described as arising from the “Fourier 

integral” 


exp 




= / dh 


exp 


if / 


exp 




(4.20) 


The above equation differs from (4.14) by factors of i due to working in Minkowski space in¬ 


stead of Euclidean space. One can experiment with how to interpret Weinberg’s statement 
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about (4.20) being a “Fourier integral”. Indeed, since e 2 is not a decaying Gaussian, one 
cannot literally perform the above integral, even its finite-dimensional analog. One could 


try to regulate (4.20) by say endowing all occurences of i in (4.20) with a small negative 
real part, i.e., replace i with i — e (this regularization has the desirable property that the 


right-hand side of (4.20) remains the exponential of a BRST exact term). This would make 


the right-hand side of (4.20) damping but then the left-hand side acquires a small positive 


real part, and so the same convergence issue reappears if one is to take the interpretation of 
integration literally. The Wick expansion however is immune to such issues of convergence, 
and is thus the actual procedure upon which the above steps relie. 

4.5. Perturbative versus constructive QFT. There is of course much more to quantum 
field theory than the perturbative approach, and so we round out our discussion with some 
brief remarks to place our perturbative analysis in a more complete context. 

To begin, we note that the value of the Wick expansion in the perturbative approach 
to a quantum field theory stems from the tacit assumption that the resulting formal series 
should provide some kind of approximation to the full theory, namely, an asymptotic series. 
As one can glean from the finite dimensional setting, this can only be true provided one can 


control the behavior of the path integral for large fields. Recall that Lemma 3.1 shows that 
the Wick expansion provides asymptotics for an integral localized via a smooth cutoff to a 


neighborhood of its critical points. Via Theorem 3.2, this boosts to an asymptotic expansion 
of the full integral given sufficient decay of the integrand at infinity. In infinite dimensions, 
establishing decay of the action in the path integral for large fields is of course much more 
nontrivial. This is quite evident in the case of nonabelian gauge theories, where gauge- 
invariance makes it difficult to establish bounds on the nonlinear terms of the Yang-Mills 
action when the gauge-field is large. 

Such large field problems pose a obstacle in relating the Wick expansion of perturbatively 
defined path integrals to the full path integral, whatever that might be. But even if one 
ignores this issue and takes the Wick expansion as the starting out, one still has the problem 
of deciding what the Wick expansion is an asymptotic series of. In finite dimensions, I{H) 
can always be made a well-defined function of h, either by taking f{x) compactly supported 
or imposing the necessary decay properties. In perturbative quantum field theory, one does 
not start off with a well-defined integral but with the formal Wick expansion defined via 
the many intermediate steps of regularization, renormalization, and possibly other tricks. 
The choice of a function whose asymptotics are governed by such a formal series is not 
unique, since e.g. functions like are asymptotically zero. Techniques for producing 

a canonical function whose asymptotics are given by a particular series (this provides a 
reverse to the diagonal arrow in Figure 1) include such methods as Borel summation. If one 
can successfully use such methods to sum a divergent asymptotic series to a functioif^ such 
a function could then be viewed as some truncated version of the path integral localized 
to fields near the Wick expansion locus. It would seem reasonable to conclude that such 
an analysis is necessary in order for one to connect purely formal perturbation theory with 
something that deserves to be called integration. See also [S] for some related perspectives 


^^This has been successfully carried out for scalar field theories in dimensions two and three [111 Chapter 
23.2]. 
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arising from “resurgence theory”, in which it is hoped that the ambigiuities that arise from 
attempting to Borel sum a perturbative series encode essential features of nonperturbative 
physics. 

We should of course note that in special cases, insight has provided suggestions of how 
to define the full path integral. We have in mind the remarkable insights of Witten, who 
(to name just two examples) used localization methods in symplectic geometry to offer an 
(alternative) construction of the Yang-Mills path integral in two dimensions |22] and insights 
from topology to define the Chern-Simons path integral for compact, oriented 3-manifolds 
m- These insights go well beyond what perturbation theory can do alone, in fact even 
bypassing what would otherwise be the analytic problem of defining measures on infinite 
dimensional space^^ or some modification thereof One of the preeminent challenges to 
placing quantum field theory on firm mathematical foundations is to make sense of the path 
integral from first principles, thereby bridging the gap between formal perturbation theory, 
which is limited to the purely formal setting of series expansions, and the few specialized 
instances in which external insights provide a serendipitous definition of the path integral. 
There are of course formulations of quantum field theory that make no use of the path 
integral, but that is outside the scope of this article. 

5. Conclusion 

We provided a basic study of the Wick expansion in finite dimensions and discussed 
its relation to integration. We then discussed how formal manipulations of perturbatively 
defined path integrals in quantum field theory are more properly understood through gen¬ 
eralization of the finite-dimensional Wick expansion and not integration. This provides 
a rigorous treatment of such formal manipulations in a conceptually sound manner. We 
concluded with a discussion of constructive quantum field theory to show that there is still 
much to be desired in making mathematical sense of the much wider features of quantum 
field theory that lie outside of perturbation theory. 
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